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The collective elementary excitations of two-dimensional magnetoexcitons in a Bose-Einstein con- 
densate (BEC) with wave vector k — were investigated in the framework of the Bogoliubov theory 
of quasiaverages. The Hamiltonian of the electrons and holes lying in the lowest Landau levels 
(LLLs) contains supplementary interactions due to virtual quantum transitions of the particles to 
the excited Landau levels (ELLs) and back. As a result, the interaction between the magnetoexci- 
tons with k — does not vanish and their BEC becomes stable. The energy spectrum contains only 
one gapless, true Nambu-Goldstone(NG) mode of the second kind with dependence u)(k) w k 2 at 
small values k describing the optical-plasmon-type oscillations. There are two exciton-type branches 
corresponding to normal and abnormal Green's functions. Both modes are gapped with roton-type 
segments at intermediary values of the wave vectors and can be named as quasi-NG modes. The 
fourth branch is the acoustical plasmon-type mode with absolute instability in the region of small 
and intermediary values of the wave vectors. All branches have a saturation-type dependencies at 
great values of the wave vectors. The number and the kind of the true NG modes is in accordance 
f^ ' with the number of the broken symmetry operators. 
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I. INTRODUCTION 

> ; 

A two-dimensional electron system in a strong perpendicular magnetic field reveals fascinating phenomena 
f* 1 **. ' such as the integer and fractional quantum Hall effects. The discovery of the fractional quantum Hall 

effect (FQHE) fundamentally changed the established concepts about charged single-particle elementary 
t^" ■ excitations in solids [1, 2]. 

^+ \ In this paper we study a coplanar electron-hole system with electrons in a conduction band and holes 

in a valence band, both of which have Landau levels in a strong perpendicular magnetic field. Earlier, 
this system was studied in a series of papers [3-9] mostly dedicated to the theory of 2D magnetoexcitons. 
This system bears some resemblance to the case of a bilayer electron system with half-filled lowest Landau 
levels in the conduction bands of each layer [10]. The coherent states of electrons in two layers happened 
to be equivalent to the BEC of the quantum Hall excitons [11] formed by electrons and holes in different 
layers. The system we are interested in has only one layer, with electrons in conduction band and holes in 
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the valence band of the same layer created by optical excitation or by p-n doping injection (both of these 
methods can be called "pumping"). In the case of a single excited layer which we consider, the density of 
excitons can be quite low, so that the electron Landau level and the separate hole Landau level are each 
only slightly occupied, and Pauli exclusion and phase space filling do not come in to play. 

II. HAMILTONIAN OF THE BOSE-EINSTEIN CONDENSATION OF MAGNETOEXCITONS. 

The effective Hamiltonian describing the interaction of electrons and holes lying on the LLLs is 

H = Hcoui + Hsuppi- (1) 
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Here Hcoui is the Hamiltonian of the Coulomb interaction of the electrons and holes lying on their LLLs: 

Hcoui = IJ2 W $ [p(Q)p(-Q) ~ Ne - fth] - Me#e - p h N h , (2) 

Q 

and Hsuppi is the supplementary indirect interactions between electrons and holes that appear due to the 
simultaneous virtual quantum transitions of two particles from the LLLs to excited Landau levels (ELLs) 
and their return back during the Coulomb scattering processes. This interaction was deduced in Ref. [9]. 
It has the form: 



ffsuppi = \b^n - ^ E v(Q)p(Q)p(-Q) - 

Q 
±-J2u(Q)D(Q)D(-Q). (3) 
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The coefficients Bi^i, V(Q) and U(Q) are proportional to the small parameter r = Ii/(huj c ) < 1, where 
/; = \/B is the ionization potential of the magnetoexciton, hw c ~ B is the cyclotron energy, and B is the 
magnetic field strength [9]. The degeneracy of the Landau level N equals S/(2irl 2 ) 1 where I is the magnetic 
length I 2 — ch/eB and S is the surface layer area. Here p(Q) are the density fluctuation operators expressed 
through the electron p e (Q) and hole f>h(Q) density operators as follows 
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p h ($) = Y,^ tl ^ Q j Qx , 

p(Q) = Pe(Q) - Ph(-Q); (4) 

D(Q) = p e (Q) + Ph(-Q); 

N e =p e (0); N h = p h (0); N = N e + N h ; 
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The density operators are integral two-particle operators. They are expressed through the single-particle 
creation and annihilation operators aLa p for electrons and b p,b p for holes. Eq is the dielectric constant of 
the background; fx e and ph are chemical potentials for electrons and holes. Coefficients V(Q), U(Q) and 
Bi-i were deduced in [9, 12]. 

As discussed in previous papers [5-9, 13], the breaking of the gauge symmetry of the Hamiltonian (1) can 
be achieved using the Kcldysh-Kozlov-Kopaev [14] method with the unitary transformation 

T(y/NZ) = cxp[^/N^(d\k) - d(k))}, (5) 

where d'(k) and d(k) are the creation and annihilation operators of the magnetoexcitons and k is the wave 
vector of the condensate. In the electron- hole representation they are [5-9]: 
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BEC of magnetoexcitons leads to the formation of a coherent macroscopic state as a ground state of the 
system with wave function 

\l> g (k)) = f{y/NZ) |0> ; a p |0) = b p |0) = 0. (7) 



Here |0) is the vacuum state for electrons and holes. Though we kept arbitrary value of k, nevertheless our 
main goal is the BEC with k — and we will consider the interval 0.5 > kl > 0. The function (7) will 
be used in the Section 3 to calculate the averages values of the type (D(Q)D(— Q)\ . The transformed 
Hamiltonian (1) looks like 
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and is succeeded, as usual, by the Bogoliubov u-v transformations of the single-particle Fermi operators 




= ua p -v(p- —)bl x _ p ; 
a p \ip g (k)) = 0; 
P P = f (v^Vel) b p fl (v^C) = ub p +v( -1 - p)a\ a 
/3 p |^(fc))=0. 



(9) 
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Here v(t) — ve lk v tl ' . The coefficients u = cos (\J2irl 2 n ex ) and v = sin (y / 2irl 2 n ex ) were determined in 
Ref.[8]. The equality v = sin(v) can be satisfied only at v < 1. 

Along with this traditional way of transforming the expressions of the starting Hamiltonian (1) and of the 
integral two-particle operators (4) and (6), we will use the method proposed by Bogoliubov in his theory 
of quasiaverages [15], remaining in the framework of the original operators. The new variant is completely 
equivalent to the previous one, and both of them can be used in different stages of the calculations. For 
example, the average values of products of two integral two-particle operators can be calculated using the 
wave function (7) and u-v transformations (9), whereas the equations of motion for the integral two-particle 
operator can be simply written in the starting representation. 

The Hamiltonian (1) with the broken gauge symmetry is written below in the lowest approximation of the 
theory of quasiavereges [15]. Side by side with the last term in (10) there are another smaller terms breaking 
the gauge symmetry. They were neglected. This approach is named as quasiaverages theory approximation 

[11]- 
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(10) 



-fjy/N (d f (k) + d(k)) 



Here the parameter fj, which determines the breaking of the gauge symmetry, depends on the chemical 
potential n and on the square root of the density similar to the case of weakly non- ideal Bose-gas considered 



by Bogoliubov [15]. In our case the density is proportional to the filling factor v 



We have: 



H = n e + n h ; jj, = n + If, N ex = v N 

fj = (E ex {k) - n)y = (E{k) - A(k) - jS)v 

E ex (k) = -I t - A(k) + E{k) 

E(k) = 2j2w QS in^ [KxQ]j2 
Q 
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A(fc) ~ r determines the shift of the ionization potential /; [9]. 



III. ENERGY SPECTRUM OF COLLECTIVE ELEMENTARY EXCITATIONS 

The equations of motion for the integral two-particle operators with wave vectors P 7^ in the special 
case of BEC of magnetoexcitons with k — arc 

i^d(P) = 
= (-fi + E(P)-A(P))d(P)- 

-2iY,W(Q>in( [Px ? ]zl2 )x 



xp(Q)d(P-Q)--x 



xJ2u(Q) 

Q 



I \P x Q]J 2 \ 
cos — - — x 



xD(Q)d(P-Q)+^ J '■ 



N 



(12) 



iH±JK?) 



= -i^W(Q)sin 



x[p(P-Q)p(Q)+p(Q)p(P-Q)}+ 

+w2^ U{Q)sm [ 2 x 



x D(P - Q)D{Q) + D(Q)D{P - Q) 



ih—D(P) 
dt v ' 



IJ2W(Q) 



/ [P x Q]J 2 \ 

sin x 



:[p(Q)D(P - Q) + D(P - Q)p(Q)] 



2N 



J2 U & 



. [PxQy?\ 

sin x 



:[D(Q)p(P-Q) + p(P-Q)D(Q)} 
d(P)-cfl(-P) 



+2ijVN 



where W(Q) = W(Q) - V(Q)/2N. 

Following the equations of motion (12) we have introduces four interconnected Green's functions 
Gij(P,t) as well as their Fourier transformations G\j(p,u>) at T = [16, 17] of the type G\j(P,t) = 
((^•(P,t);B(P,0))), where Ai(P,i) - d(P,t),A 2 (P,t) = d<<(-P,t),A 3 (P,t) = p{P,t) / y/N , M{P,t) = 
D(P, t)/vN and an arbitrary operator B(P, 0) because its choice does not influence on the energy spectrum 
of the system. These Green's function can be named as one-operator Green's functions. At the right hand 
sides of the corresponding equations of motion there is a second generation of the two-operator Green's 
functions. A second generation of equations of motion derived for them containing in their right hand sides 



the three-operator Green's functions. Following the procedure proposed by Zubarev [17] the truncation of 
the chains of equations of motion was made and the three-operator Green's functions were presented as a 
products of one-operator Green's functions Gij(P,u>) multiplied by the averages of the type (D(P)D(-P)). 
The average values were calculated using the ground state wave function (7) and the u-v transformations 
(9). For example, it was obtained 

(D(P)D(-P)) = 4u 2 v 2 N; (p(P)p(-P)) = 0. (13) 

The Zubarev procedure is equivalent to a perturbation theory with a small parameter of the type v 2 (l — v 2 ). 
The closed system of Dyson equations has the form 

4 

^G li (P,w)V ij (P,w) = C lj ; 3 -1,2,3,4. (14) 

There are 16 different components of the self-energy parts Y,jk{P,oj) forming a 4 x 4 matrix. Four diagonal 
self-energy parts are: 
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They contain real and imaginary parts as follows Ejj(P, w) = <Tij{P, lo) + iT(P,u>). The cumbersome dis- 
persion equation for electron-hole collective excitations is expressed in general form by the determinant 
equation 
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Sy(Aw) =0. (16) 



Due to the structure of the self-energy parts, it separates into two independent equations. One of them 
concerns only the optical plasmon branch, corresponding to oscillations of the difference of electron and hole 
densities and has a simple form: 

E 33 (P,w)=0. (17) 

It does not include at all the chemical potential p, and the quasiaverage constant i). But the average values 
depends also on the condensate wave function (7). 



The solution of the equation (17) is 
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The right hand side of this expression at small values of P has a dependence \P\ and tends to saturate 
at large values of P. The optical plasmon branch Hcoop(P) has a quadratic dispersion law in the long 
wavelength limit and saturation dependence in the range of short wavelengths. Its concentration dependence 
is of the type -\/v 2 (l — v 2 ) what coincides with the concentration dependencies for 3D and 2D plasmons. 
The obtained dispersion law is represented in the Fig.l. 

The second equation contains the self-energy parts £n, £22, £44, £14, £<4i, £24 an< ^ ^42) which include 
the both parameters \x and 77 and has the form 

S 11 (P;w)S 22 (P;w)i;44(P;u;) - 

-E 4 i(P;^)S22(P;w)S 14 ( J P;^)- (19) 

-£42(P;w)£ii(P;w)£ 24 (P*;u;) = 0. 

The solutions of the equation (19) describe the exciton energy and quasienergy branches arising due to the 
normal and abnormal Green's functions as well as the acoustical plasmon branch. The ideal magnetoexciton 
gas can exist only in the case v 2 = 0, with an infinitesimal number of excitons, but without plasma at all. 
In this case the plasmon frequency is zero whereas the exciton of one exciton is possibly. 

In the zero order approximation on the concentration, when v = v 2 = the real parts au(P,uj) look as: 

o-ii (P,w) = -<7 22 (-P, -w) = 

= hw + jL-E(p V ) + A{P)- 

<j 33 (p',u J )=<j 44 (p',lo) = 0- (20) 

^=0;,u + A(0) = 

In this limit the energy needed to excite the magnetoexciton with k — in the state with P / equals 
to E(P), whereas the plasma oscillations do not exist at all. In this approximation the exciton branches 
of the spectrum are gapless and true NG modes. In the first order approximation on the small parameter 
v 2 (l — v 2 ) we have a non-ideal Bose gas and the self-energy parts (15) contain terms linear and quadratic 
in the interaction constant U(P). The last terms have the denominators with unknown frequency what 
increase the number of the solutions. They also can be taken into account also by the iteration method. 
The first step in this direction gives the following parts <Tij(P,u>) 



a 41 (P,u;) = -fj + U(P)^; (21) 



a n (P, u) = -a 22 (-P, -«) = hui + p. - E(P) + A(P); 

, <d(0)> 

^ { P,u)=fj-U(-P)W§; 
<ju(P,u) = —27); <j 2 a{P,u) = 2fj;a iA (P,uj) = hu>; 

Now the exciton branches of the spectrum transform from the true into the quasi-NG modes with the gaps 
in the point P = as follows 



hu = ± a 



,- B (P, + A,»f^,-^|S) ,22) 



In this approximation the acoustical plasmon branch continues to vanish. Taking into account the quadratic 
terms on the interaction constant U(P) we have obtained the corrections to the exciton branches of the 
spectrum and the dispersion law for the acoustical plasmon mode. 
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FIG. 1. Three branches of the collective elementary excitations: the exciton-type quasi-NG mode with a gap in 
the point pi — (dash-dotted line); the second-type NG mode describing the optical plasmons(solid line) and the 
first-type NG mode with absolute instability (dotted line) describing the acoustical plasmons (dashed line). 

The acoustical plasmon branch corresponding to oscillations of the total particle density has a dispersion 
law completely different from the optical plasmon oscillations. It has an absolute instability beginning with 
small values of wave vector going on up to the considerable value PI ~ 2. In this range of wave vectors, 
the optical plasmons have energies which do not exceed the activation energy t/(0). The supplementary 
Hamiltonian gives rise to general attraction in Hartree approximation, characterized by the coefficients 
U(0). It plays the role of activation energy if one should like to overpass this attraction. The ground state 
of the system is unstable as regards the generation of the acoustical plasmons. It means that the system 
becomes a localized generator of the growing acoustical waves without their traveling through the system 
as in the case of convective instability. 

In the case of 2D magnetoexcitons in the BEC state with small wave vector kl < 0.5 described by the 
Hamiltonian (10), the both initial continuous symmetries are lost. It happened due to the presence of the 
term rj(cft( k ) + d{ k )) in the frame of the Bogoliubov theory of the quasiaverages. Nevertheless the energy 
of the ground state as well as the self-energy parts Sy(P, ui) were calculated only in the simplest case of 

the condensate wave vector k =0. These expressions can be relevant also for infinitesimal values of the 
modulus |k| but with a well defined direction. In this case the symmetry of the ground state will be higher 
than that of the Hamiltonian (10), what coincides with the situation described by Georgi and Pais [18]. It is 
one possible explanation of the quasi-NG modes appearance in the case of exciton branches of the spectrum. 
Another possible mechanism of the gapped modes appearance is the existence of the local gauge symmetry, 
the breaking of which leads to the Higgs effect [19]. The interaction of the electrons with the attached 
vortices gives rise to a gapped energy spectrum of the collective elementary excitations as was established 
in Ref.[20, 21]. The number of the NG modes in the system with many broken continuous symmetries was 
determined by the Nielsen and Chadha [22] theorem. It states that the number of the first-type NG modes 
Ni being accounted once and the number of the second type NG modes TV// being accounted twice equals or 
prevails the number of broken generators Nbg- It looks as follows Nj + 27V// >z Nbg- In our case the optical 
plasmon branch has the properties of the second-type NG modes. We have Nj = 0; Njj = 1 and Nbg = 2. 
It leads to the equality 27V// = Nbg- The three branches of the energy spectrum are represented together in 
the Fig.l. One of them is a second-type Nambu-Goldstone(NG) mode describing the optical plasmon- type 
excitations, the second branch is the first-type NG mode with absolute instability describing the acoustical- 
type excitations and the third branch is the quasi-NG mode describing the exciton-type collective elementary 
excitations of the system. It is a gapless true NG mode in zero order approximation on the small parameter 
v 2 (l — v 2 ) and become gapped in the first order approximation due to the fact that the symmetry of the 
ground state with small condensate wave vector k is greater than the symmetry of the Hamiltonian. 



IV. CONCLUSIONS 



The energy spectrum of the collective elementary excitations of a 2D electron- hole system in a strong per- 
pendicular magnetic field in the state of Bose-Einstein condensation with wave vector k = was investigated 



in the framework of the Bogoliubov theory of quasiaverages. The starting Hamiltonian describing the e-h 
system contains not only the Coulomb interaction between the particles lying on the lowest Landau levels, 
but also a supplementary interaction due to their virtual quantum transitions from the LLLs to the excited 
Landau levels and back. This supplementary interaction generates, after the averaging on the ground state 
BCS-type wave function, direct Hartree-type terms with an attractive character, exchange Fock-type terms 
giving rise to repulsion, and similar terms arising after the Bogoliubov u — v transformation. The interplay 
of these three parameters gives rise to the resulting nonzero interaction between the magnetoexcitons with 
wave vector k = and to stability of their BEC as regards the collops. 

The separated electrons and holes remaining on their Landau orbits can take part in the formation of 
magnetoexcitons as well as in collective plasma oscillations. Such possibilities were not taken into consid- 
eration in the theory of structureless bosons or in the case of Wannier-Mott excitons with a rigid relative 
electron-hole motion structure without the possibility of the intra-series excitations. 

The energy spectrum of the collective elementary excitations consists of four branches. Two of them 
are excitonic-type branches (energy and quasienergy branches). The other two branches are the optical 
and acoustical plasmon branches. We can say that results obtained in our system are similar to what 
was obtained in system of BEC of the quantum Hall excitons [11]. In these both models there is only one 
gapless Nambu-Goldstone mode between four branches of the energy spectrum. In both models the exciton 
branches of the spectrum are not gapless and differ from the NG modes. In our case the exciton energy 
and quasienergy branches corresponding to normal and abnormal Green's functions have a gaps in the 
point P = 0, a roton-type segments in the region of intermediary wave vectors PI ~ 1 and saturation-type 
behaviors at great values of PI. 

We should like to thank the referee, whose remarks stimulated to essentially improve our paper. 
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